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Abstract 

We show how Wick polynomials of random variables can be defined combinatorially as the 
unique choice which removes all “internal contractions” from the related cumulant expansions, 
also in a non-Gaussian case. We discuss how an expansion in terms of the Wick polynomials 
can be used for derivation of a hierarchy of equations for the time-evolution of cumulants. 

These methods are then applied to simplify the formal derivation of the Boltzmann-Peierls 
equation in the kinetic scaling limit of the discrete nonlinear Schrodinger equation (DNLS) 
with suitable random initial data. We also present a reformulation of the standard perturba¬ 
tion expansion using cumulants which could simplify the problem of a rigorous derivation of 
the Boltzmann-Peierls equation by separating the analysis of the solutions to the Boltzmann- 
Peierls equation from the analysis of the corrections. This latter scheme is general and not 
tied to the DNLS evolution equations. 
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1 Introduction 

Wick polynomials, also called Wick products, arose first in quantum field theory as a way of 
regularizing products of field operators. The principal goal there was to replace monomial products 
by polynomials with state dependent coefficients, chosen so as to remove singular terms appearing 
in the associated perturbation expansion. 

The procedure can also be applied in more general probabilistic settings. The following def¬ 
inition is given in Wikipedia [1] and in the Encyclopedia of Mathematics [2]. Consider n (real) 
random variables yj, j = 1,2,..., n, on some probability space (fl, y) and denote expectation 
over the probability measure y by (•). The Wick polynomial with powers kj > 0, j = 1, 2,..., n, 
are then defined recursively in the total degree ki + k 2 + ■ • ■ + kn by the following conditions: 

1. li ki = k 2 = ■ ■ ■ = kn = 0, set = 1. 

2. If the total degree is greater than zero, require that {'-yi^y^^ ■ • ■yn''-) = 0. 

3. For all j, require that the (algebraic) derivatives of the Wick polynomials satisfy 

dy, ■ ■ • yf • • • •' • yf~" • • • . (1.1) 


These conditions have a unique solution for which ly^^y^^ ■ ■ -yn"' is a polynomial of total degree 
ki + k 2 + ■ • ■ + kn in the variables yj. (The uniqueness is algebraic, not only almost everywhere as 
random variables. That is, the conditions fix all coefficients of the polynomials. This can be seen 
by induction in the order |/|: the requirement in item 3 fixes all new coefficients apart from the 
constant, which is then fixed by the vanishing of the expectation value in item 2.) The coefficients 
are polynomials of expectations of the random variables yj, and hence depend on the measure fx. 
The first order polynomial is obtained by simply centering the variable, :yi: = yi — {yi), but already 
at second order more complex structures appear, :yiy 2 . = yiy 2 — {yi)y 2 — {y 2 )yi — {yiy 2 ) + ‘^{yi){y 2 )- 
If the random variables have joint exponential moments, i.e., if there is /3 > 0 such that 
(e^^jl^^l) < oo, the Wick polynomials can also be obtained by differentiating a fairly simple 
generating function. It can then be defined for A G R", such that |Aj| < (3 for all j, by 


G-w{X',yi, ..., 2/n) 


exp (Efci 

(exp(EEiA.2/.))’ 


and then for all kj > 0, j = 1,2,... ,n. 


■■yi^y2^ 


•: = O': 




Gw(A;yi,..., yn) 


A =0 


( 1 . 2 ) 


(1.3) 


The generating function Gw (A; yi,..., yn) is also called “Wick exponential” and often denoted by 
“:exp (Er=i Xiyi):”. For a derivation and basic properties of such Wick polynomials, see [3]. 

The Wick polynomials become particularly simple to use if the joint measure of y is Gaussian. 
Defining the covariance matrix by Cj>j := Cov{yji,yj), a Gaussian measure has (exp(A.7/)) = 
gA.(y>+A.CA/ 2 _ Therefore, the generating function of the Wick polynomials then reads simply 
Gw(A; y) = exp[A • (y — (y)) — A • GA/2], and the resulting Wick polynomials are closely related to 
Hermite polynomials. This is the setting encountered in the original problem of renormalization of 
quantum field theories (the “unperturbed measures” concern free fields and hence are Gaussian). 

In the Gaussian case, one can also identify the Wick polynomials as arising from an orthogonal- 
ization procedure. Wiener chaos expansion and Malliavin calculus used for stochastic differential 
equations can be viewed as applications of such orthogonal projection techniques. [4] 

In the non-Gaussian case, there are far fewer examples of applications of Wick polynomial 
techniques. The computations become then more involved. For instance, there is no explicit 
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formula for the generating function unless the inverse of the moment generating function happens 
to be known explicitly. In addition, then the polynomials typically no longer form an orthogonal 
set in 

The goals of the present contribution are two-fold. In the first part, we show that Wick poly¬ 
nomials have a natural combinatorial definition, closely connected to cumulants and the related 
cluster expansions of correlation functions. We also rederive their main properties without resort¬ 
ing to the generating function, hence without assuming Gaussianity or the existence of exponential 
moments. 

In the second part, we show how Wick polynomial expansions may be used in the analysis of 
stochastic processes. In particular, the goal there is to apply the expansion to study the time- 
evolution of the cumulants, i.e., of the connected correlation functions, of the process. We will 
explain there why often it is cumulants, and not moments, which should be used as dynamical 
variables. For simplicity, we consider here only processes whose dynamics are deterministic and 
given by a differential equation, such as Hamiltonian evolution in classical particle systems. The 
randomness enters via the initial state. However, generalization to Markovian stochastic dynam¬ 
ics should be straightforward, for instance, if the generator of the process maps polynomials to 
polynomials. 

In the general setup, the best one can hope for are recursion relations leading to an infinite 
hierarchy of equations connecting the evolution of the cumulants. We explain in Section 4 what 
immediate constructions are available for hierarchical study of the evolution of cumulants and 
Wick polynomials. 

We give more explicit applications in section 5 where we study the evolution on a lattice of 
particles following the discrete nonlinear Schrodinger (DNLS) equation with random initial data. 
In particular, our goal is to show how the Wick polynomial expansion of the dynamics greatly 
simplifies the (still only formal) derivation of the related Boltzmann-Peierls equation. This case is 
one of the few examples of nonlinear Hamiltonian evolution where a rigorous analysis of the related 
perturbation expansion has been possible so far: in [7] it was proven that if the initial measure is 
a stationary Gibbs measure, then the time-correlations of the field follow an evolution equation 
derived using a perturbation expansion analogous to the one needed for the Boltzmann-Peierls 
equation. 

An ultimate goal of the present reformulation of the evolution problem would be to complete 
the rigorous derivation, and hence give a region of validity, of the Boltzmann-Peierls equation. We 
show how the Wick polynomial expansion could help in this goal by separating the problem of 
solving the effect of the Boltzmann-Peierls evolution from the estimation of the corrections arising 
from the wave nature of the microscopic evolution, such as constructive interference. For the 
DNLS evolution the Wick polynomial expansion coincides with what was called “pair truncation” 
in [7]. In fact, the present work arose from an attempt to generalize this construction to other 
polynomial potentials, which we later realized to coincide with Wick polynomial expansions. 

It should already be apparent from the above example that in order to use the Wick polynomials 
some care is needed in the choice of notations to avoid being overcome by lengthy formulae and 
intractable combinatorial estimates. We begin by explaining our choices in detail in Section 2. 
The first part containing the combinatorial definition and properties of Wick polynomials is given 
in Section 3. The second part discussing the use of Wick polynomial expansions for the study of 
evolution of cumulants begins in Section 4. We conclude it with the specific application to DNLS 
dynamics in Section 5. Some comments and possible further directions are discussed in Section 6. 
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2 Setup and notations 

We consider a collection yj, j G J where J is some fixed nonempty index set, of real or complex 
random variables on some probability space If yj are complex, we assume that the 

collection is closed under complex conjugation, i.e., that to every j there is f G J such that 

y^’ = y*3- 

Expectation over the probability measure y, will be denoted by E or (•). In case the underlying 
measure needs to be identified, we denote the expectation by E^ or (•)^. We use sequences of 
indices, I = (ii, 12 ,..., G J", to label monomials of the above random variables, with the 
following shorthand notation 


= y^ly^2■■■y^n^Y[y^k■ (2-i) 

k^l 

We also set y® := 1 if / is the empty sequence. Since all yj commute with each other, we have 
y^ = y^ for any two sequences I, I' which differ by a permutation. 

We will need to operate not only with such sequences but also with their subsequences and 
“partitions”. This will be done by choosing a distinct label for each member of the sequence 
and collecting these into a set. How the labelling is done is not important, as long as one takes 
care when combining two “labelled” sets. We rely here on the following standard conventions: 
any sequence (ik) can be uniquely identified with the function k iu which itself is uniquely 
determined by its graph, the subset {{k, ik) \ k = 1,2,... ,n} of N x J. We consider subsequences 
to be subsets of the graph of the sequence. Partitions of the sequence then correspond to partitions 
of its graph which can be understood as partitions into nonempty subsequences. 

Mathematically, this leads to the following structure. Finite (sub)sequences of indices are now 
uniquely labelled by the collection J^, which consists of those finite subsets H C N x J with the 
property that if {n,j), {n',j') G A and {n,j) ^ {n',j') then n ^ n'. We also allow the sequence 
to be empty which is identified with 0 G J^. For nonempty sets, the natural number in the first 
component serves as a distinct label for each member in A. In addition, we can use the order 
of the natural numbers to collapse any H G back to a sequence H in J: Given H G with 
n > 0 elements, there is a unique bijection g : {1, 2,..., n} A such that its first component is 
increasing, g{k)i < g{k')i for all k < k'. Using this g, we define Ak := g{k )2 G J for fc = 1, 2,..., n. 

To each finite sequence I = (ik) of n elements in J, we assign I := {{k, ik)\k = 1,2,... ,n} 
as the set of labels. Obviously, then I and any of its subsets belong to J^. The following list 
summarizes some basic notations and definitions which will be used later without further remark. 

1. If / is a sequence, and a set is needed by the notation, the set is chosen to be /. For instance, 
the notation “H C /” means A C I. 

2. The notation ViE) denotes the collection of partitions of the a set U G J^. If / is a sequence, 
V{I) ■.= V{I). 

3. If H G and it is used in a place of a sequence, the formula always refers to the collapsed 
sequence A obtained via the increasing bijection g above. For instance, then y^ := y^ = 
rifc^i J/g(fc) 2 - (Note that if / is a sequence, then y^ = y^ in agreement with (2.1).) 

4. If H G J^, we denote the corresponding sequence of random variables by yA '■= (VAk)^^!- 

5. If TO G N and H G J^, the notation refers to a set where any element with label to is 

cancelled, i.e., := {(k,ik) € A\k ^ to}. Note that it is possible that = A. If / is 

a sequence and to < |/|, /I™! corresponds to a sequence which is obtained by removing the 
TO:th member from I. 

6. If I, I' are two sequences, they can be merged into a new sequence (U,..., i\i\,i'i,..., i\ji\) 

which we denote by / + 11 A,B G J^, we take A + B := A + B. For a merged sequence. 
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the notation “/ C I + /'” always refers to the collection of the labels of the first |/| members 
and analogously C / + refers to the collection of the last |/'| members. The merge 
operation is clearly associative, and we hence drop parentheses when it is applied iteratively; 
for instance, / + /' + I" is a sequence of length |/| + |/'| + \I"\. 

7. To avoid separate treatment of expressions involving empty sets and conditions, we employ 
here the following standard conventions: if the condition P is false, we define 

^(•••):=0, and set also P(0) := {0}. (2.2) 

p p 

Similarly to the moments, to any / € we denote the corresponding cumulant by one of the 
following alternative notations 

K[yi] = ■ ■ ■ ^,y^„] = ■■■ ,yij- (2.3) 

The corresponding Wick polynomial is denoted by 

■y^ly^2■■■yin■ = ^V■ = ^V■^i ■ (2-4) 

Note that this notation is slightly formal, since the result is not a function only of the power 
y^ but depends on all subpowers, y^, A C I, a.s well. It also requires that one carefully defines 
which random variables are being “Wick contracted”. We will use parentheses for this purpose, if 
necessary. For instance, “:{y^):” means y^ — E[y^] which usually differs from :y^:. 

As an application of the above definitions, let us point out that the earlier defining Wick 
polynomial condition (1.1) is equivalent to the requirement that for every nonempty sequence I 
and any j € J we should have 


Id 

dypV-= 

fc=i 

Here, and in the following, 1 denotes the generic characteristic function: 1(P) = 1 if the condition 
P is true, and otherwise l(T’) = 0. 

We recall that, if the random variables yj, j = 1, 2,..., n, have joint exponential moments, then 
moments, cumulants and Wick polynomials can be generated by differentiation of their respective 
generating functions which are 

G„,(A) := E[e^-], g,{X) := \nG^{X) and Gw(A; j/) := ^^ = . (2.6) 

Here X • x := X]r=i clarity we have denoted the integrated random variable 

by “x” instead of “y”) and the “generation” happens by evaluation of the /:th derivative at zero, 
i.e.. 


E[/] = 9{Gm(0), K[yi] = digdO) and :/: = (9 {Gw( 0; y), (2.7) 

where “i9{” is a shorthand notation for . 

As a side remark, let us also recall that it is possible to replace the above definitions of 
generating functions by parametrizations which do not require the existence of any moments and 
hence work for arbitrary Borel probability measures g. If all yj are real, then replacing the 
exponential e'^'^ by e''^'^ yields an L^{y) function for all A € R". If y G C", the same is achieved 
by using e^-^ ■v+>'-y )/2 ^nd A G C": in this case, differentiation with respect to ReA^- generates 
“iReyj” and with respect to Im Aj generates “ilmy^”. Naturally, without absolute integrability of 
the moments it is not guaranteed that any of the derivative exist. However, it might nevertheless be 
useful to inspect the time evolution of the generating function, in particular, if the time evolution 
is regularizing and improves the integrability of the moments. 
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3 Combinatorial definition and properties of the Wick poly¬ 
nomials 

Let us first recall the “moments-to-cumulants” formula which holds for any I G as long as all 
moments y^, A C I, belong to 

E[/] = XI n , (3.1) 

77^V{I) AGtt 

where V{I) denotes the collection of partitions of the set I. (Or to be precise, of I. Here 
it is important to assign a distinct label to each random variable in the power to get the 
combinatorics correctly.) For a partition tt S 'P{T), let us call the subsets A G tt clusters or blocks. 
Let us also recall that the cumulants are multilinear, i.e., they are separately linear in each of 
the variables yj, j G I. These results are discussed, for instance, in [5] and [6] and also briefly in 
Appendix A here. 

Let us point out that by the conventions adopted here, (3.1) is indeed valid also for the empty 
sequence / = 0. Then the sum over partitions is not empty since it contains tt = 0. However, the 
corresponding product is empty since there is no A with A G tt. Therefore, the right hand side of 
(3.1) evaluates to one which agrees with our definition of E)?/®]. 

We next show that it is possible to choose a subset of the indices and remove all its “internal 
clusters” from the moments-to-cumulants formula by replacing the corresponding power with a 
polynomial of the same order. This will be achieved by using the following recursive definition. 

Definition 3.1 Suppose that Iq G is such that E[|y^|] < oo for all I C Iq- We define polynomi¬ 
als W[y^] := /o’’ C /o inductively in |/| using the following rule: set W[y®] := 1, 

and for / yf 0 use 

E Eiy''] ■ (3.2) 

The definition makes sense since the W-terms on the right hand side all have an order lower than 
|/|. It also implies that indeed each W[y'^] is a polynomial of order |/|, with only the term y^ 
being of the highest order. It is also straightforward to prove by induction that the coefficients 
ceW] can be chosen so that they only depend on with A C I. Therefore, the definition 

of W[y^] is independent of Iq in the following sense: if IqjIq G are such that < oo for 

every I C Iq and I C Iq, then for all / C Iq D Iq we have W[y^; Iq] = /g]. In Appendix A 

we explain how cumulants can also be defined via a similar recursive construction. 

The following theorem shows that these polynomials indeed have the promised truncated 
moments-to-cumulants expansion. We also see that the polynomials are essentially uniquely de¬ 
fined by this property. What is perhaps surprising is that the coefficients of the polynomial can be 
chosen depending only on the moments of its constituent random variables. This implies that the 
same polynomial can be used for many different probability distributions, as long as the marginal 
distributions for the constituent random variables are the same. 


Theorem 3.2 Assume that the measure y has all moments of order N, i.e., suppose that E[|2/'^|] < 
oo for all I G with |/| < N. Use Definition 3.1 to define W[y^] for every such I. 

Then replacing y^ by W[y'^] removes all terms with clusters internal to I: the following trun¬ 
cated moments-to-cumulants formula holds for every I' G with |/'| -I- |/| < 




W[y ]y 


7rG'P(/+/') 


1(A n yf 0 for all A G e) E,[yA\ . 


AGtt 


(3.3) 


In particular, E[W[?/^]] = 0 if I $. 

In addition, if I G with |/| < N/2 and W' is a polynomial of order at most |/| such that 
(3.3) holds for all I' with |/'| < N — \I\, then W is y-almost surely egual to yV[y^]. 
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Corollary 3.3 Assume that E[| 2 /^|] < oo for all I G ^. Then W[y^] are ^-almost surely unique 
polynomials of order |/| such that (3.3) holds for every I' GJ^. 


Proof: We make an induction in |/|. By (3.1), the claim is true for |/| = 0 since then J = 0 and 
thus W[y^] = 1. 

Assume then that I ^ % and that the claim is true for sets of size less than |/|. Consider an 
arbitrary /' G such that |/'| + |/| < N. For E C I, denote := (/ + 1') \ E. Given a partition 
TT of / + /', we can define 

TTi := {Ai G TT I Ai n /' 7 ^ 0} 


and 


TTo := 7r\7ri . 


Then E := Utto C I and ttq G V{E), tti G V{E‘^) (also whenever E or E‘^ happens to be 
empty). Once tt is fixed, the decomposition tt = ttq U tti is unique and we thus find that 1 = 

= ^0 UttOHaig^i O/' 7 ^ 0). Using this in the standard 
moments-to-cumulants formula shows that 




= E 


.,I+T 


Y. n 


= Y Y Y Y 1 ( 7 ’' = TTo UtTi) Y[ l^lVAo] n {K[yA^]^Al n r ^ 9)) 

7 re(P( 7 +/^) EdI 7 TodV{E) 7 TidV{E^) Aq^ttq AiGtti 


= E E EH n MvAMAini'^m 

EGI 7ToGV{E) 7TiGV{E‘^) Aq^TTq A\d7T\ 


= ^ 1(A n/'7^ 0 VA G tt) 

irGP(/+/') A^-k 


+ E E n '^[22^0] E 1 (^0/'7^ 0 VA G TTi) Y [ K [ yA ,] 

0^E(Zl ttoGE{E) Aq^ttq TTiGEiE*^) Ai^tti 




1(A n I' 7 ^ 0 VA G tt) kIpa] 

AGtt 


E 

<D^E(ZI 


E[y^]E\w[y^\V 


where in the last step we used the moments-to-cumulants formula and the induction hypothesis 
(note that E’^ collapses to the sequence {I\E) +1'). Hence, by the definition (3.2) equation (3.3) 
holds for this I. This completes the induction step and shows that (3.3) is valid for all I, /' with 
|/| -I- |/'| < N. If / 7 ^ 0 and /' = 0, we have / -I- /' 7 ^ 0 so that for any tt G P(/ -I- 1') there is some 
A G TT and then obviously Afil' = %. Thus (3.3) implies that E[W[y^]] = 0 for / 7 ^ 0. 

To prove uniqueness, suppose that / G with |/| < N/2 and W' is a polynomial of order at 
most |/| such that (3.3) holds for all I' with |/'| < N — |/|. Then Pj := W' —W’[y^] is a polynomial 
of order at most |/| and E[Piy^ ] = 0 for all I' with |/'| < N/2. Since the collection of random 
variables is assumed to be closed under complex conjugation, this implies that also E[P/(y^ )*] = 0 
whenever |/'| < N/2. Thus we can take a linear combination of such equations and conclude that 
E[|P/p] = 0. This implies that Pj = 0 almost surely, i.e., that W' = W[y^] almost surely. 

This concludes the proof of the Theorem. The Corollary is then an immediate consequence. □ 


In fact, the polynomials given by Definition 3.1 are equal to the standard Wick polynomials. 


Proposition 3.4 Suppose that Iq G is such that E[|y^|] < 00 for all I C Iq- Then W[y^] = :y^: 
for every I C Iq- 

Proof: If / = 0, we have W[y^] = 1 = :y^:, and else by Theorem 3.2 we have E[W[y^]] = 0. 
Therefore, to prove W[y^] = :y^: it suffices to check that (2.5) holds when the Wick polynomials 
are replaced by W-polynomials. We do this by induction over |/|. Firstly, if / = 0, we have 
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dyjy^W] = 0, as required. Assume then that / ^ 0 and that the claim is true for sets of size less 
than |J|. For every j € J and nonempty E C I, the induction assumption implies that 

FI 

dy,W[y^\^] = Y.Mzk = j)H{k,ik) ^ (3.4) 

fe=i 


(The second characteristic function allows keeping the labeling inherited from / by adding zero 
terms into the sum for the “missing” labels.) Since W[y'^] satishes (3.2), we thus find that its 
algebraic derivatives satisfy an equality 




\i\ 

= ^M^k = j) 


j(k) 


E[y®]l((fc,4) ^ A)>V[/''‘’\-®] 

0^ECI 


FI 

= Y 



^ E[y^]W[y^'\^] 


/'=/('=) 


1^1 


(3.5) 


where in the last equality we have applied (3.2). This proves that also satisfies (2.5) 

and hence completes the induction step. Therefore, the polynomials W[y^] satisfy the defining 
properties of Wick polynomials and thus W[y^] = :y^:. □ 


3.1 Basic properties of the Wick polynomials 

In this section, we assume that there is Iq G ^ is such that IE[|y^|] < oo for all I C Iq- This 
guarantees the existence of the Wick polynomials W[y^] for all / C Iq, and allows using the 
results from the previous section. In particular, by Theorem 3.4 these are equal to the standard 
Wick polynomials and from now on we will use the standard notation :y^: for them. 

The next Proposition collects some of the most important properties of Wick polynomials. 

Proposition 3.5 The following statements hold for any I C Iq: 

1 . 


2 . 

3. 


^Y -y^' = Y -y^' Y 11 • 

UCI UCI TveV{I\U) AeTr 


(3.6) 


Wick polynomials are permutation invariant: if I' is a permutation of I, then :y^ : = :y^:. 


V-- = Yy'' Y 

Ucl 7TeP{I\U) AGtt 

4- If I' := denotes the sequence obtained by cancelling the first element of I, 


(3.7) 


V- = yii -V- - Y '^[2/v] = yii -V- - Y ^[2f(n)+c/] ■ (3-8) 

(1M)&VCI UCI' 

Proof: Item 1: The first equality in (3.6) follows directly from the definition (3.2) since then 

y^ = Yny^\ ^ . 

Ed UCI 
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(3.9) 



The second equality follows then by using the moments-to-cumulants expansion. 

Item 2: The permutation invariance of the Wick polynomials follows using straightforward 
induction in the definition (3.2) since the random variables commute and hence the powers are 
always permutation invariant. 

Items 3 and 4: Let us first define W[y^] by setting it equal to the right hand side of (3.7) for any 
I C lo- If / = 0, wehave >V[ 2 /^] = 1 = :y^:. Suppose / ^ 0. Since X(^eP( 0 )(-l)'’"' OyieTr 
the definition yields a polynomial of order |/| in y. Our goal is to prove that 

_ \i\ _ 

dy,my'] = E =J)W [/'“’], (3.10) 

k^l 

=yiiyv[y^"^] - K[yi]- E . (3.11) 

(i,*i)eyc/ 

Then the claim = :y^: follows by straightforward induction in |/|: Case |/| = 0 was proven 

above. Suppose / 7 ^ 0 and that W[y^ ] = :y^ : whenever |/'| < |/|. Then the induction assumption 
and Theorem 3.2 can be used to evaluate the expectation of the right hand side of (3.11), implying 
E[>V[y^]] = 0. By (3.10), the polynomial W[y^] also satisfies the third defining condition of the 
Wick polynomials, equation (2.5). Hence, W[y^] = :y^: which completes the induction step. Then 
(3.11) implies the first identity in (3.8) and the second identity is found by a relabeling of the 
summation variable. Hence, also item 4 follows. 

To prove (3.10), consider some / 7 ^ 0. In the definition of W[y^], we can express the derivatives 
of y’^, U C I, as in (3.4). This shows that 

/ 

dyjW[y^] = '^^Hik=j)mk,ik)&U)y^^'‘^ ^ (-l)l’"l K[yA] 

Ucl k=l TrGV{I\U) AGtt 

\I\ 

= Ei(*fc = j) E (-I)'’"'n 

fe=l 7reP(/<'“)\V) 

\I\ _ 

= ^l(zfc=j)W[/'‘E (3.12) 

fe=i 

Therefore, (3.10) holds. 

To prove (3.11), denote x := (l,ii) and /' := I^^\ We first split the definition into two parts 
as follows: 

^[ 2 ^^]=E ^ E (- 1 )'’"'n +E ^ E (- 1 )'’"'n ^[ 2 /a] ■ 

Ucl TT<SV(I\U) A^tx ucl TXCV(I\U) ACtx 

(3.13) 

Following a reasoning similar to (3.12), we find that the first term in the sum on the right hand 
side is equal to yi^ W [j/'^ ]. The second term is equal to 

Y,l{x^U)y^ ^ (-l)H(' ll ^[y^] 

ucl TrC'P{I\U) ^X^AC-IX 

= - E 1 ( 2 ; ^ u)y^ E '^[2 /v] E (-I)'’"' n ^[ 2 /^] 

ucl xCVCl\U TXCV((I\U)\V) ACtt 

= - E "^[22^] E 2/^ E (“I)'’"' n 

XCVCI UCI\V TxCV(iI\V)\U) ACtx 

= - E i'^{yv]'^[y^'^^] ■ (3-14) 

xCVCl 
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Therefore, (3.13) implies that also (3.11) holds. This completes the proof of the Theorem. □ 
Example 3.6 Written is terms of cumulants, the Wick polynomials of lowest order are 
:y: = y- K{y), 

■yiy2- = yiy2 - K{yi,y2) - K{yi)y2 - K{y2)yi + K{yi)K{y2), 

■yiy 2 y 3 ' = 2/11/22/3 - k (2/ i ,2/2,2/3) + k (2/ i ,2/2) k (2/3) + K{yi,y 3 )K{y 2 ) + (2/2,2/3) k (2/1) 
- k (2/ i ) k (2/2) k (2/3) - k (2/ i .222)2/3 - k (2/ i ,2/3)222 - ^(2/2,2/3)221 
+ k ( i / i ) k (2/2)2/3 + n{yi)n{y 3 )y 2 + K{y 2 )K{y 3 )yi 
- k (2/ i )2/2223 - K{y2)yiy3 - n{y3)yiy2 ■ 

Proposition 3.7 The Wick polynomials are multilinear, i.e., if a,/S are constants such that yj = 
ayi + Pyr for some j, i, i' G J, then, whenever I and k are such that ik = j, we have 

:y^\ = a :y^^ *+*^*): +/3 -.y^^ . 

Proof: The claim follows using multilinearity of cumulants in the representation formula (3.7). □ 

The following result extends the earlier theorem and shows that multiple application of Wick 
polynomial replacements continues to simplify the moments-to-cumulants formula by removing all 
terms with any internal clusters. 


Proposition 3.8 Assume that the measure n has all moments of order N, i.e., suppose that 
E[|22^|] < 00 for all I € with |/| < N. Suppose L > 1 is given and consider a collection of L + 1 
index sequences J' ,Ji G , £ = 1,... ,L, such that |d'| \Ji\ < N. Then for I + 

(with the implicit identification of Ji and J' with the set of its labels in I) we have 


E 


Y[y^:y^' 

i=i 


E n {K[yA]l{A (f Ji \/£)) . 


(3.15) 


Proof: We proceed via a double induction: the first induction is over L and the second induction 
is over \Jl\. The case L = 1 follows directly from Theorem 3.2. Now we assume as induction 
hypothesis of the first level induction that L >2 and 


L-l 

n 

£=1 




E n (k[7/a]1(^ t Ji'J£< L)) . 

7reP(Efci' P+J') 


(3.16) 


Then we consider the second induction over \Jl\ ='■ m- For m = 0, we have '.y^^: = 1 and thus 
then the induction hypothesis (3.16) directly implies (3.15). As induction step of the second level 
hypothesis we take that, for fixed L, equation (3.15) holds for all \.Jl\ < m. Then, if \Jl \ = m, we 
can use (3.2), (3.16) and the second level induction hypothesis to justify the following argument 


10 




analogous to the one used in the proof of Theorem 3.2 


E =E ^ E[y^] E 

^1=1 -I L^=i -I IH^EcJl '-^=1 

E n {K[yA]HA (^Ji'il< L)) 

7reP(Efci^ Jt+JL+J') 

" E EH E n 

ID^EgJl -iro&ViE) Aeno Je + (JL\E) + J') A'Gtti ^ 

X 1(A (;tJeW< L)1{A' Jl\ E)^ 

= E n iK[yA]l{A^JeW)) 

7rG'P(/) A^tt 

+ ^2 1(3A € TT s.t. A c Jl) (K[y^]l(A JiWi < L)) 

7rG'P(/) 

- Y. n'J{A&E\AcJL]^^)\{{E[yA]1{AtJi^^<L)) 

7reP(/) .4671- 

= E n {K[yA\l{A t Jl V£)) . 

•7reP(/) 


This completes the induction step and hence also the proof. 


4 Cumulants and Wick polynomials as dynamical variables 

To study the state of a random system, it is often better to use the cumulants rather than the 
moments of the random variables. For instance, if j/, z are independent random variables we 
have E[7 /"z’"] = E[j/"]E[z’"], which is typically nonzero, whereas the corresponding cumulant is 
zero whenever both n, m ^ 0. Hence, for systems where two regions become “asymptotically 
independent” (for instance, for sufficiently mixing stochastic processes), it is the cumulants, not 
moments, which will vanish in the limit. 

To have a concrete example, let us consider a random lattice field yx, x a Z"^, which is spatially 
sufficiently strongly mixing. Then, for instance, K{yo,yi) —> 0 if the distance of the index set 
/ C from the origin becomes unbounded. Often in the applications the mixing is so strong that 
the cumulants, in this case also called connected correlation functions, become ‘Ti-clustering”: 
for any order n one requires that sup 2 ,g 2 d \Eiyx,yi)\ < oo. Naturally, such a property 

should then not be expected to hold for standard correlation functions E[j/^], apart from some 
degenerate examples. 

In addition to being mixing, the random fields found in the applications are often translation 
invariant. This means in particular that all moments E[j/^] remain invariant if every index in I is 
translated by a fixed amount, i.e., = E)?/^] for every a; G Z"^ if we set I{x)e := ii — x. If the 

system is both £i-clustering and translation invariant, the cumulants of the Fourier transformed 
field yk := k indexed by the d-torus T'^, satisfy 

n 

K[yikuk2,...,k^)] = s(^ke^Fniki,k2,...,kn), ( 4 . 1 ) 

where “6” denotes the Dirac delta distribution and the arithmetic on is defined via periodic 
identification. Here F„ denotes the Fourier transform of F„(V) := l{Xi = 0)K{yx), X G (Z'^)”, 
and for £i-clustering measures Fn is a uniformly bounded continuous function of k. Therefore, 
although the cumulants are singular, their singularity structure is simple, entirely encoded in the 


II 



i5-multiplier. In contrast, by the moments-to-cumulants formula, then for / := (1,2, ...,n) and 
any k € (T'^)” 

E n (4.2) 

T7eV{I)A&TT ^ JGA ' 

which has ever more complicated singularity structure as the order of the moment is increased. 
(The above discussion can be made mathematically rigorous by replacing the infinite lattice by a 
periodic d-dimensional lattice. See [7] for details.) 

Hence, for stochastic processes which lead towards a state which is mixing and translation 
invariant, it seems better to focus on the time-evolution of cumulants instead of the corresponding 
moments. However, it is not immediately clear how to avoid resorting to the moments as a middle 
step. It turns out that using Wick polynomials instead of “bare” monomials to define the time- 
evolution helps in achieving this goal. Recall that any monomial term can always be expanded 
in terms of Wick polynomials using (3.6), albeit with state-dependent coefficients. 

To have a concrete example how this could work in practice, we consider in the following the 
case of deterministic evolution with random initial data. Explicitly, we assume that the system 
is described by random variables j/j(t), where j belongs to a fixed (finite) index set J, and t > 0 
denotes time. The initial values yj{0) are random with a joint distribution /xq, and for each 
realization of y(0) we assume that the values at later times t > 0 are determined from the solution 
to the differential equation 




(4.3) 


where the functions Mj{t) are “interaction amplitudes” from the /:th Wick polynomial of y(t) 
to yj{t). For each j G J, the set collects those / G which have a nonzero amplitude, i.e., 
Mj (t) 7 ^ 0 for some t > 0. For simplicity, we assume here that is finite and that the amplitudes 
Mj{t) are some fixed functions of time. (They typically might depend on the cumulants of y{t), 
but this is not relevant for the discussion below: it suffices that they are not random variables.) 

We present a concrete example of such a dynamical system in Appendix B where we show 
how the evolution of N classical particles interacting via a polynomial interaction potential can 
be described by a system of this type assuming a known random distribution of initial positions 
and momenta. Another explicit example is given in Section 5. 

The usefulness of representing the dynamics in the form (4.3) becomes apparent when looking 
at the evolution of cumulants. To avoid technical complications, let us suppose that the joint 
exponential moments of y{t) exist and are continuously differentiable and uniformly bounded 
functions of t. This will simplify the discussion since it allows using the generating functions 
defined earlier in (2.6). With the shorthand notation X-y := ^jVj, the cumulant generating 

function of y{t), t fixed, is gt{X) = ln(e'^'^(*)), and the Wick polynomial generating function is 
G{X;y{t)) = The time-evolution of the cumulant generating function is 

connected to the Wick polynomial generating function by 

dtgtiX) = (e^'*'(‘))-i(A • dty{t)e^-y‘^*'>) = (A • dty{t) GiX;ym . (4.4) 

Therefore, for any /' 7 ^ 0 (using the slightly symbolic notations “/ \ i” and “j/j” instead of 
and yi^, when i = {k, ik) G I) 

dtK[y{t)i'] = 9ydtgtiX)\^^^ = '^{^ty^{t) :y{ty '\":). (4.5) 

Hence, if the evolution satisfies (4.3), we obtain 

dtK[yit)r] = E E ■y^f'^"-) ■ (4-6) 

ie/' le^i 
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In this case, determining the evolution of expectation values of all multiplications of two Wick 
products, where both Ii,l2 are non-empty, would also yield a solution to the 

evolution of cumulants. 

We can now obtain a closed cumulant evolution hierarchy using (4.6) and Theorems 3.2 and 
3.8. First, note that for any / 0 and j € J we have 

{:y{ty: -yjit):) = {:y{ty:yj{t)) = K[y(t)/+(y)], (4.7) 

since in this case there is exactly one non-internal cluster, the entire set I + (j). In addition, if 
/ = 0, we clearly have {:y{t)^: :y{tY :) = 1(7' = 0). Therefore, the first two cumulants satisfy, for 
arbitrary j, f G J, 

dtn[yYt)\ = 1(0 e yj)M^{t ), (4.8) 

dtK[y{t)u,j')] = M/(t)K[y(t),+(,,)] + ^ M/,(t)At[j/(t)^+(,)]. (4.9) 

For higher order cumulants, with |/'| > 3, the equation typically becomes nonlinear; we then have 

dtK[y{t)i'] = ^ ■^i)A7‘f^^(t)K[2/(t)o)+(/'\i)] +J2 ^(y)^-yyy'■vyy''^"-) ■ 

iei'jeJ ier iej^i,\i\>2 

(4.10) 

We have separated here the terms with |/| = 1 to show how they operate linearly on the cumulants 
of order |7'| (note that K[j/(i)(j)-i-( 7 '\i)] = i^[y{t)i"] fo'' the sequence I" which is obtained from I' 
by replacing i with j). In the final sum, both |7| and \I' \ i\ are greater than one, so it has a 
cumulant expansion 

{■.y{ty::y{ty'\^:)= ^ H ^ 7 ^ 0, A n (/'\ *) ^ 0)) , (4.11) 

T^eV{i+(r\i)) agtt 

i.e., all clusters have to contain at least one element from both sequences. In particular, it cannot 
contain any singlets, i.e., it does not depend on any of K[yj{t)], j G J. Let us also point out that 
since in these terms |7| -|- |7' \i| > |7'|, any linear term is necessarily of higher order. In particular, 
this means that lower order cumulants can appear only in nonlinear combinations. 

Instead of studying the full cumulant hierarchy, one can also use evolution estimates for the 
Wick polynomials. The situation often encountered in the applications is that the properties of the 
initial measure are fairly well known, whereas very little a priori control exists for the time-evolved 
measure. In such a case, one can use the above result and obtain a perturbation expansion by 
applying the fundamental theorem of calculus. With the shorthand notation y := y(0) we have 

K[y{t)i'] = K[yi']+ [ Ml{s){:y{sy::y{sy'\^:) 

“'0 ier lejfi 

= i^[yi'] + ^^{-y^--y^''^"-) j dsMl{s) 

ier leJ^i 

+ J2J2 f ds'ds'{-yis'y---y{s'y''^'') [ dsMlis), (4.12) 

iev leJ’i ds' 

where we have applied Fubini’s theorem to the final integral. This type of expansion could be 
helpful if the coefficients (isM/(s) behave better than M/(t), such as in the presence of fast 
oscillations. Further iterations of this procedure, using either the above cumulant hierarchy or 
any of the Wick polynomial hierarchies below, would then yield an expansion of K[y(t)ii\ in terms 
of the expectations at time t = 0 and the time dependent amplitudes This is particularly 

useful if all it) are small, since each iteration adds one more such factor. 
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Let us conclude this section by deriving recursion formulae for the products of Wick polyno¬ 
mials. As mentioned earlier, these could then be used instead of the direct cumulant hierarchy to 
study the time-evolution of the cumulants. For this, it would suffice to study 
appearing in (4.6), but typically the products of two terms do not satisfy a closed evolution equa¬ 
tion and a full hierarchy will be needed. Let us begin with the evolution equation for For 

any deterministic evolution process, we can obtain a fairly compact evolution equation by treating 
the time-derivative dtUj as a new random variable: 

dt'.yitY: = didtG{X;y{t))\^^^ = ^ . (4.13) 

iei 

The form is analogous to the standard Leibniz rule. For products of Wick polynomials, we thus 
have 

n n 

9t n ■yYY''- = H ■yYY'^'- ■ (4-i4) 

k—1 k—1 i£lk k'^k 

Assuming (4.3) and using multilinearity, we then obtain the following equation involving “nested 
Wick products”: 

- n 1 n r 

dtE ■{■yYY-)yYY'’''"-Y[-yitY’‘'- ■ (4-i5) 

-k—1 k—li^IklG<^i ^ k'^k 

The formula (4.15) is appealing in its simplicity but it does not directly lead to closed hierarchy 
of equations. This can be achieved by expanding the nested product in terms of cumulants and 
Wick products. For this, let us note that by (3.8) and the observation made after (4.7), we have 
for any /' 

■■{dtyr{t))y{tY’■■ = dtyi{t) :y{tY'■■ - ^ ^^ty^{t):y{tY’■]:y{tY''^'^■ ■ (4.16) 

ucr 

Therefore, whenever (4.3) holds, we find that 

n n f n+1 

9tE[n:y(i)''“:] =EE E ^^^(t) E[:2/(t)^'=\b 

fc=l k=l ielk ^ k' = l;k'^k 

r " nr n 1 

- E n ■yitY’’'-\^[-y{tY"^^--y{tf-\\ (4.17) 

k' = l\k' ^k 

This forms a closed hierarchy of evolution equations for the collection of all expectation values of 
the type E OLi ■vitY'"- ■ 

A second alternative for the hierarchy follows from the observation that if y{t) and z{t) are 
two processes which start with independent, identically distributed initial data, then at any later 
moment they are also independent and identically distributed and hence 

dtG{X; y{t)) = G(A; y{t)) (A ■ dty{t) - E^[A • dtz{t) G(A; z{t))]) (4.18) 

= E;^ [G(A; y(t))G(A; z{t))X ■ {dty{t) - dtz{t))] 

where in the second equality we have used Ez[G(A; z{t))] = 1. Consider then the product measure 
for the processes y, z and let G' denote the corresponding Wick polynomial generating func¬ 
tion. Since by Fubini’s theorem then E„,^[e^4!/(t)-HAt))] = E^[e^ 2^(*)]2, now G(A; y(t))G(A; z(t)) = 
G'{X;y{t) + z{t)). Hence, for dynamics satisfying (4.3) 


dtEy 



n 

EE E M/(t)E,.,[(:y(t)^:-:z(t)^:) -.(yit) + z{t)Y^\^-. J] ■ 

k—1 iGlk I^<^i k'^k 

(4.19) 
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Let us point out that the earlier expression in (4.17) follows from the above one if we expand the 
power {y{t) + and then use the fact that the joint measure is a product measure. The 

formula does not yet yield a closed hierarchy but the following generalization does so: if Zk,i{t) 
are processes such that their joint initial distribution is given, then 




n'(E 

fc=i I 


-EE EE^/w 

k=liGlkI&yi i 


Ik\i 


k'^k t 


Zk',l' 


(4.20) 


where refers to a measure where Zk/ for each £ has been independently duplicated in the 
z'k ^-process. 

A possible benefit of this formulation could be when Zk,i{t) have mean zero and are independent 
for all Then the central limit theorem governs the behavior of Zk,i when there are many 
terms in the sum. Therefore, it could be of help in controlling the otherwise difficult case where 
one has performed many iterations starting from (4.19). 


5 Kinetic theory of the discrete NLS equation revisited 

In this section we apply the previous Wick polynomial techniques to the discrete nonlinear 
Schrodinger equation. This example is chosen since it has a particular simple, but nontrivial, 
Wick expansion of the evolution equation. In addition, we can then rely on the rigorous results 
and known properties from an earlier work on the model [7]. We focus on the kinetic theory of 
the model on the “kinetic” time-scale which is 0{\~Y if the nonlinear coupling A is taken to zero. 

We begin by going through the derivation of the Boltzmann-Peierls equation in the spatially 
homogeneous case, and show how the task is simplified by using the Wick expanded dynamics 
and the cumulant hierarchy, as explained in Section 4. We only consider terms which would be 
present in the simplified case of Gaussian initial data. We give an example in Appendix C which 
highlights the mechanism leading to suppression of the additional effects of non-Gaussian initial 
data in the kinetic scaling limit. 

For this particular setup, it is easy to find dynamical variables whose evolution equation does 
not have a linear part. This is an important simplification since it negates a term which becomes 
rapidly oscillating on the kinetic time-scale, having a divergent frequency 0(A“^) in the kinetic 
scaling limit. The effect becomes apparent when looking at field time-correlations instead of the 
evolution of equal time cumulants. We discuss the issue in more detail in Section 6.1. 

In [7], the initial data is taken to be given by a thermal Gibbs measure which is stationary 
both in space and in time. We do not assume the initial data to be time-stationary here, but the 
computations in this section require space translation invariance. The spatially inhomogeneous 
case is technically substantially more complicated, and we discuss it only briefly in Section 6.1. 

The results in this section are derived in the spirit of standard perturbation theory and focus 
solely on evolution on short kinetic time scales, t = rA“^ with 0 < t 1. It is however possible to 
apply the cumulant hierarchy differently, leading to equations which do not require taking a scaling 
limit. We conclude the study of the DNLS model in Section 6.2 by proposing a reformulation of the 
problem which leads to Boltzmann type evolution equations which could be accurate also for times 
longer than 0(A“^). The discussion is not completely mathematically rigorous, but we propose 
a scheme which could be used to this end under some natural conditions about the time-evolved 
state. 

The discrete NLS equation on the lattice deals with functions T/f : IR x > C which satisfy 
idtYY-x) = ^ a(x - y)V’t(y) + A|'!/'t(x)|Vt(x). (5.1) 

yGZ'^ 
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Here the function a : ^ IR is called the hopping amplitude and we assume that it is symmetric, 

a(x) = a(—x), and exponentially decreasing. The parameter A > 0 is considered to be small, and 
in the kinetic scaling limit we take A —>■ 0 and t = with r > 0 fixed. The initial field 'ipo 
is assumed to be random, bounded on finite subsets of the lattice, and to have an £i-clustering 
distribution. We aim at controlling the moments of the random variables '0t(x) and which 

we label using + 1 ) := V’i(x) and — 1 ) := V’t(x)*. 

Since we do not assume that tpo is t^ 2 -summable, even the (almost sure) existence of solutions 
to (5.1) becomes an issue. To our knowledge, it has not been proven for the above setup, and 
most likely, some additional assumptions about the increase of the values of the initial field at 
infinity are needed for proper existence theory. However, these problems can be easily avoided 
by replacing the infinite lattice by a finite lattice with periodic boundary conditions: see [7] 
for details. This would merely result in replacing the lattice, the Fourier space and transform, 
and the associated ^-functions by their finite lattice counterparts. Since even then the final limits 
cannot be rigorously controlled, we opt here for some additional formality in the discussion, but 
with less complicated formulae to deal with. 

For technical simplicity, here we also only consider initial data which are “gauge invariant”: 
we will always suppose V'o(x) has the same distribution as e‘^^/>o(x) for any 9 G [0, 27r]. In fact, 
this transformation commutes with the time evolution, i.e., if the initial field is changed from 
tpo to e'^' 00 ) the time-evolved field will change from ipt to In particular, also the field ipt 

will then be gauge invariant. The main reason for insisting on this assumption is that it will 
automatically force many cumulants to be zero and hence simplify the combinatorics associated 
with the cumulant hierarchy. Gauge invariance implies that a moment is zero unless it has the 
same number of ■0 and ip* factors, even when the fields are evaluated at different times. Hence, it 
implies that every odd moment of the fields is zero and hence also every odd cumulant. Similarly, 
we see that even cumulants are also zero if they concern a different number of ip* and p) variables. 

For instance, we find using (3.7) that for any gauge invariant state and any Oj := iptji'^jjO'j), 
J = 1,2,3, 

:aia 2 a 3 : = 010203 — £[ 0102)03 — £[ 0103)02 — £[ 0203)01 . (5.2) 

This is the definition of the “pairing truncation operation” V given in Lemma 3.2 of [7). Applying 
the truncation operation in the evolution equation was one of the key changes to the standard 
perturbation theory which allowed the rigorous analysis in [7]. With the benefit of hindsight, we 
can now identify it as a Wick contraction of the random variables. 

Under the above assumptions and using (3.6), we find that (5.1) is equivalent to the following 
Wick contracted evolution equation 

idtp’ti^) = X! «(x-y):0t(y):+2A£['0t(x)*'0t(x)]:'0t(x):-|-A:0>t(x)*0t(x)V>t(x): . (5.3) 

Hence, the random variables ipti^, cr), a = ±1, satisfy an evolution equation of a form required in 
the previous section, in (4.3), 

i(jdt'tpti^,cr) = ^ a(x-y):V’t(y, cr):-l-Ai?t(x):'0t(x, cr):-|-A:'0t(x,-l)'0t(x, ( 7 )' 0 t(x, I):, (5.4) 

where we have defined Rt{x) := 2£[|'0t(x))^) > 0, which is also equal to 2K('0t(x, —l),'0t(x, I)). 

5.1 Translation invariant initial measures 

The evolution problem simplifies significantly, if we assume that the initial data is not only gauge, 
but also translation invariant. Since also spatial translations commute with the time evolution, 
we can use the earlier results for the cumulants of Fourier transforms of the random field. In 
particular, then for any t > 0 and x, y G Z"^, 

£[0>t(x)*'0t(y)) = £[0>t(x - yY'ptiO )]. (5.5) 
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This implies that i?t(x) = i?t(0) =: Rt for all t, and therefore the evolution equation (5.4) for 
translation and gauge invariant initial data can be written as 

iadtipt{x,cr) = at^(x-y):'!/'t(y,cr):+A:'!/'t(x,-l)7/>t(x,CT)V’t(x, 1):, (5.6) 

where a^(x) := a(x) + Al(x — 0)Rt- 

Using multilinearity of the Wick polynomials, we thus find the following evolution equation for 
the Fourier transformed fields ^t(k, ct) := e“‘^’^*' ’^'0t(x, cr), 

dtiptiKcr) = -iCTWj^(k) :V't(k, cr): 

-iaX [ (iki(ik2dk3(5(k - ki - k2 - ks) :'0t(ki,-I)^t(k2,cr)^t(k3,1):, (5.7) 

J(T‘^)3 

where 


a;^(k) := S^(k) = S(k) + \Rt. (5.8) 

For later use, let us point out that the definitions of the random fields imply the following rule for 
complex conjugation of the Fourier transformed fields: ^/’t(k, ct)* = —a). In addition, the 

assumed symmetry of a implies the symmetry a;^(—k) = a;^(k). 

We recall that in the present translation invariant setting, the n:th cumulants satisfy for k G 

(T‘^)", a e {-1,1}” 


K(^t(ki,CTi),... ,^t(k„,CT„)) = s(^Y^^)^riik,a;t), (5.9) 

where F„(x,a;t) := l(xi = 0)«:(^/>t(xi, ui),..., ^/>t(x„, cr„)) is identically zero unless '^gcri = 
0. We are now mainly interested in the evolution of the lowest nonzero cumulants, i.e., of 
F 2 (x, (—1,1)). We denote its Fourier transform by W; more precisely, we set 

W,\k) ^ = E ■ (5-10) 

It follows that F 2 ((ki,k 2 ), (-1, l);t) = lF/'(k 2 ) and F 2 ((ki, k 2 ), (1,-1); t) = lFt^(-k 2 ). There¬ 
fore, we have the following general rule for second order cumulants: 

«:(^t(ki,cri),^t(k2,cr2)) = (5(ki -|- k2)l(CTi + a2 = 0)lU/((T2k2). (5-11) 

Therefore, to study the evolution of all second moments in this systems, it suffices to study 
the function W^. In particular, clearly Rt = 2 Jj^dk.W^(k). Note that by translation invariance, 
is always real valued. 

5.1.1 Heuristic derivation of the Boltzmann-Peierls equation 

After these preliminaries, we are ready for an application of the cumulant hierarchy to study the 
Fourier transformed fields. Our first goal is to justify the Boltzmann-Peierls equation which has 
been conjectured before, based on perturbation expansions in [7]. The conjecture says that in the 
kinetic scaling limit the function W should converge: it is assumed that there exists a limit 

IIA(k) := lim IFA_ 2 (k). (5.12) 

In addition, the analysis of the perturbation series suggests that the limit satisfies the following 
homogeneous Boltzmann-Peierls equation: 

dtWtik) = C{Wt{-)){k), (5.13) 
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with the collision operator 


C(l^(-))(k) =47r 


/(Jd)3 


(iki(ik2(ik3(5(k + ki — k2 — k3)(5(w + wi — a;2 — W3) 


X [Vh(ki)Vh(k2)ty(k3) + Wik)Wik2)W{ks) - W{k)W{ki)W(ks) - W{k)W{ki)Wik2)] 

(5.14) 

where uji := a{ki) and w := S(k). 

In fact, a lucky accident hides the fact that our present random fields are actually ill suited for 
taking of the scaling limit: it is clear from the linear part in (5.7) that they are highly oscillatory, 
and only observables where these oscillations cancel out, can be hoped to have a (nonzero) limiting 
value in the kinetic scaling limit. Fortunately, there is a simple “renormalization” which cancels 
these fast oscillations. If we define a new random field by the formula 


at(k, cr) = iptik, a) exp ( icr 


ds ujt 


(k) 


(5.15) 


then it clearly satisfies an equation without a linear term. Explicitly, then 


dtat{ki,a) = -iaX / 


dk2dk3dk46{ki — k2 — k3 — k4) 


X -l)at(k 3 , a)atik4, 1): . 


(5.16) 


Note that due to the alternating signs, the time dependent terms cancel each other out in the 
oscillatory phase term inside the integral. In fact, the same happens in the second order cumulants, 
as can be checked by using (5.11) and the symmetry of we then find that 

K(at(ki,(Ti),at(k2,CT2)) = (5(ki +k2)l(cri + a2 =0)W^{a2k2). (5-17) 

It is clear that multiplication with a nonrandom term as in (5.15) does not spoil the gauge invari¬ 
ance of the held so we can rely on it also when working with the cumulants of the a-helds. 

We can now study the evolution of W^ik) by employing the expansion given in (4.12) to the 
cumulant K{at{k', a'), at(k, cr)). We then hnd using any a' = —a that 

5(k' -h k){Wt^{ak) - W^{ak)) 

= — iAcr f ds f dkidk2dk3(5(k — ki — k 2 — k 3 )e*^^'^‘^''"“^“'^‘^^“‘^^^«:[(ao)/] 

Jo J{jd)3 

- iXa' f ds f dkidk2dk3Sik' - ki - k 2 - k3)e‘"('"'‘^'+“i-‘"'‘^"-“==)Av[(ao)7-] 

Jo J{j<i)3 

/ ds [ dkidk2dk3(5(k-ki-k 2 -k3)e*"('^‘^+‘^i-'"“'=-‘^^) 

Jo Jij-^r 

X f ds' f dkidk^dkeSiki - k4 - kj - :af/: ] 

Jo J{j<i)3 

- AV' V / ds f dkidk 2 dk 3 Sik' - ki - k 2 - k3)e*"('"'‘^'+“i-‘"'‘^=-“3) 

Jo d(T<i)3 

X f ds' f dk4dk5dk6<5(k7-k 4 -k 5 , (5.18) 
Jo J{j<i)3 

where 


/ = ((ki,-l),(k2,a),(k3,l),(k',a')) , 

/' = ((ki,-l),(k 2 ,a'),(k 3 ,l),(k,a)) , 
Je = ((k 4 , - 1 ), (kg, 0 - 7 ), (kg, 1)) . 
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(5.19) 

(5.20) 

(5.21) 



Following the standard perturbation recipe, we next apply the cumulant hierarchy to the terms 
depending on as' in ( 5 . 18 ). This resnlts in a sum of two terms: one, in which every a^' has been 
replaced by oq, plus a “correction” which we denote by d(k' + k)7?.3(crk, t). Further iterations of 
the perturbation expansion and a careful study of the oscillations of the term by term expansion 
as in [ 7 ] leads us to the conjecture that TZsCk, tA~^) should converge in the kinetic scaling limit, as 
A —> 0 , at least for sufficiently nondegenerate dispersion relations and for large enough dimension 
d. In addition, the analysis indicates that the limit value is O(t^), which is negligible compared to 
the contribution from the other terms following from ( 5 . 18 ). However, the term by term analysis 
does not suffice to actually prove the claim since the method which was used to rigorously control 
the convergence of the perturbation expansion in [ 7 ] was based on time stationarity of the initial 
state. This assumption cannot be made here since we are interested in nontrivial time evolution 
effects. Instead of going into the details of the above argument, we discuss a less technically 
involved motivation for the claim in Section 6 . 2 . 

Thus now we need to evaluate expectations of the form £[1010203: :040506:] where each at 
stands for one of the field variables. The cumulant expansion in Theorem 3.8 and the vanishing 
of the third order cumulants imply 

E[:ai0203: :040506:] = «;(oi, 04)«;(o2, 05)«:(o3, 05) + /«(oi, 04)/«(o2, 06)«:(o3, 05) 

+ «;(oi,a5)«;(o2, 04)«:(a3, 05) + «:(oi, a5)s;(o2, 06)«:(a3, 04) 

+ «;(oi,a6)a;(o2, 04)^(03,05) + «:(oi, a6)«;(o2, 05)«:(o3, 04) 

+ «;(oi,02,03,04,05,05) + “9 X K2K4' ( 5 . 22 ) 


where the last contribution denotes a sum of the nine terms consisting of a product of a second 
order cumulant and a fourth order cumulant. Naturally, also some of the above terms can be zero 
because of the gauge invariance constraints. 

To better work with the expressions arising from ( 5 . 18 ), let us next introduce a few shorthand 
notations. We denote 


/ 


(I(k 


IF(k) 

:= <(k). 


dki2...„ 

:= / d]iidk2 

••dk„ , 


1 

1 

II 

-kfc). 


:= w(ki)-1-a;(k2) 

1 

CO 

1 


'll 

1 

to 

- w(k3) -|-u;(k) 


( 5 . 23 ) 


We also choose ct = I, ct' = — I and we will only consider the pairing contractions (i.e., the Gaussian 
contractions) in the expansion ( 5 . 22 ). In fact, all terms arising from the non-pairing contractions 
are typically negligible in the kinetic scaling limit of the present type. As an example, in Appendix 
C we show how the first order terms in ( 5 . 18 ) vanish in the kinetic limit by assuming sufficient 
regularity of the dispersion relation w and the £i-clustering property of the fourth order cumulants. 
As explained in [ 7 ], the contribution from the non-pairing terms in ( 5 . 22 ) can be controlled by 
similar techniques but we will skip this more involved analysis here. 

Hence, after integrating out the variables k^, i = 4 , 5 ,6, the fourth term in ( 5 . 18 ) gives 


2A^d(k-f k') 


dki23^(k - ki23)IF(k)IF(k2)IF(k3) 







dki 23 (I(k - ki 23 )IF(k)IF(-ki)IF(k 3 ) 
dki23(I(k - ki23)IF(k)IF(-ki)IF(k2) 


dki 23 (I(k - ki 23 )IF(-ki)IF(k 2 )IF(k 3 ) 


ds 


J 

0 

ds 


J 

0 

■t 

ds 



/o ^0 


+ NPC 


( 5 . 24 ) 
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where NPC stands for ’’non-pairing contraction terms”. We proceed in the same way for the fifth 
term in (5.18) yielding 


2A2(5(kk') f dki23<5(k + ki23)lP(k)W(-ki)W(-k2) f ds f 

iJ Jo Jo 

- /dki23<5(k + ki 23 )lP(k)W(-ki)W(k 3 ) f ds f 

J Jo Jo 

- f dki23<5(k + ki 23 )lP(k)W(-k 2 )W(k 3 ) f ds f 

J Jo Jo 

+ f dki23<5(k + ki 23 )lP(-ki)W(-k 2 )W(k 3 ) f ds f 

J Jo Jo 


-hNPC. (5.25) 


By changing integration variables so that ki —> —ki in (5.24) and ki —> —k 3 , k 2 —s- —k 2 , k 3 —ki 
in (5.25), we obtain 

Wt^(k) - W(k) - (7e3(k,f) +NPC) 

= 2X‘^ f dki235(k-fi ki - k 2 - k 3 ) / ds f 
J{j'i)3 Jo Jo 

X [W(k)W(k2)W(k3) - W(k)W(ki)W(k3) - W(k)W(ki)W(k2) -f W(ki)W(k2)W(k3)] 

-f2A2 f dki23(5(kki - k 2 - k 3 ) f ds f 

J Jo Jo 

X [W{k)W{k2)W{k3) - I4^(k)V^(ki)H^(k3) - W{k)W{ki)W{k2) + V^(ki)H^(k2)V^(k3)] 

rt 


= 2A2 f dki23<5(k ki - k 2 - k 3 ) f ds f dr 
J Jo J-s 


X [W(ki)W(k2)W(k3) -h W(k)W(k2)W(k3) - W(k)W(ki)W(k3) - W(k)W(ki)W(k2)]. 

(5.26) 


Note that for any fl G IR we have f* ds dr e*”^ = dr e‘’'^(A^t — A^|r|). By setting 

t = tA“^ and taking A 0, this expression formally converges to r f^dre'^^ = T27rS(ft). 
Therefore, doing this in (5.26) yields the conjecture that 


Wr(k) — W(k) — 0(t^) = rdTT J dki23d(k -|- ki — k2 — k3)d(u; -|- wi — u;2 — oJo) 

X [W(ki)W(k2)W(k3) -h W(k)W(k2)W(k3) - W(k)W(ki)W(k3) - W(k)W(ki)W(k2)]. 

(5.27) 

Since here W(k) = ITo(k), if we divide the left hand side by r and then take r ^ 0, it converges 
to drWri^) at T = 0. Dividing the right hand side of (5.27) by r yields C(Wo(-))(k), as defined 
in (5.14). Therefore, the Boltzmann-Peierls equation (5.13) should hold at r = 0. Assuming that 
the state of the original system remains so regular that the estimates leading to the conjecture 
continue to hold, we thus find that the Boltzmann-Peierls equation should be valid for the limit 
of also at later times t, as was claimed in the beginning of the subsection. 


5.1.2 Decay of field time-correlations 

As a second example of how the standard perturbation expansion works for the cumulants, we 
consider a kinetic scaling limit of time-correlations. In particular, our goal is to show how the 
main results proven in [7] relate to the present cumulant hierarchy expansions. 

The notation “dj” was used in [7] to define the finite periodic lattice analogue of the present 
Ot-field. (One can compare the definition of “dt” in (3.9) and its evolution equation in (3.10) in [7] 
to those given for a* here.) Translated to the present infinite lattice setup, the main theorem of [7] 
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(Theorem 2.4) leads to the following conjecture about the decay of time correlations of at'- start the 
system from an £i-clustering equilibrium Gibbs state. Then there is a continuous function A^(k) 
such that E[ao(k', —l)at(k, 1)] = (j(k' + k)A^(k). The conjecture is that the kinetic scaling limit 
of exists and its decay is governed by the “loss term” of the Boltzmann-Peierls equation (5.13) 
evaluated at the corresponding limit equilibrium covariance function lT®'^^(k) = /3“^/(w(k) — /x) 
where /? > 0 and /x € IR are parameters determined by the equilibrium state. (Such functions 
are indeed stationary solutions of (5.13).) More precisely, [7, Theorem 2.4] is consistent with the 
conjecture that 

lim (k) = , (5.28) 


where 


r(lb(-))(k) = -2 / dr f dkidk2dk35(k + ki - k 2 - 

Jo J(j'i)3 


/(T'i) 

X [lT(k2)lT(k3) - lT(ki)lT(k3) - lT(ki)lT(k2)] 


(5.29) 


Instead of assuming that the system starts from an equilibrium state, let us consider more 
general states which we assume to be gauge and translation invariant and £i-clustering. We can 
immediately use the results derived in Section 4 if we consider the “oq” term to be a new field 
which has trivial time evolution with zero amplitudes, i.e., the corresponding J^j-set is empty. The 
net effect of this change is that more than half of the terms analyzed in the previous section will 
be absent. For instance, the expansion of the time correlation using (4.13) reads 

s;[ao(k', -l),at(k, 1)] 

= K[ao(k', -1), ao(k, 1)] - J ds j dki 23 (I(k - ki23)e'®(‘^+“i““^““^^K[(ao)/'+(k'-i)] 

- A2 ^ / ds [ dkidk 2 dk 3 (j(k - ki - k 2 - k3)e'"(‘^+“i 

X I ds' J dk4dk5dk65(k^-k4-k5-k6)e‘"'('"'^“^+“"-'"^‘^"-‘^'‘)E[:af/::ao(k',-l)af;''’:] 

(5.30) 

where /' = ((ki,-1), (k 2 ,1), (k 3 ,1)) and Ji = ((k4,-1), (ks, cr^), (ke, 1)). 

As in Section 5.1.1, we now assume that only pairings contribute in the kinetic scaling limit. 

When applying (5.22) to expand E[:ay : :ao(k', —l)a^/ :], we note that every pairing term results 

in a product containing a factor A^, (—k') and a product of two bF^Vtcrms. The rest of the structure 
is identical to the one considered earlier, some of the terms are merely missing now. We then use 
the perturbation expansion to the product once more. This produces a term where s' is set to 
0, and a remainder which we assume to be negligible as before. The rest of the computation is 
essentially the same as in Section 5.1.1, yielding 

A:(l;,^_ 2 (k) — AQ(k) — (terms higher order in A or r) 

= 2X^1 dki23<5(k + ki - k 2 - ka) / ds [ 

d(T‘i)3 Jo Jo 

X [Ao^(k)Wo^(k2)Wo^k3) - Ao\k)Wo^(ki)Wo^k3) - A^(k)Wo^ki)Wo^k2)]. (5.31) 

Hence, if we divide the equation by r and then take A —0, followed by r —> 0, we find that 
Ar{)d) := limA-i-o A:^^_ 2 (k) should satisfy at r = 0 

drArik) = -A,(k)r(W,(-))(k), (5.32) 

where T has been defined in (5.29). As before, the conjecture is that this equation continues to 
hold for other values r > 0, as well. 
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Once Wr is given, equation (5.32) is straightforward to solve. Since ^o(k) = Wo(k), the 
solution reads 


yl^(k) = . (5.33) 

If the system is started in an equilibrium state with Wq = 11^®'^*, we have II4 = for all s. 

Thus (5.33) implies (5.28) in this special case. 


6 Discussion about further applications 

6.1 Limitations of the direct renormalization procedure: inhomoge¬ 
neous DNLS 

The field renormalization used with the translational invariant data greatly simplihed the evolution 
equation by removing the linear term. The renormalization procedure, given in (5.15), was a simple 
multiplication by a time and fc-dependent function and the time-dependent first order terms had no 
effect in the interaction term. Unfortunately, this case is atypical: most commonly, the necessary 
renormalization is not a multiplication operator and the first order terms will also affect the 
oscillatory phase terms arising from the harmonic evolution. In fact, this happens also for the 
DNLS model as soon as we drop the requirement that the initial data is translation invariant. To 
explain the changes needed in the renormalization procedure, we discuss in this subsection the 
DNLS model with inhomogeneous initial data in some more detail. 

Before considering the inhomogeneous case, let us begin with an example which emphasizes the 
importance of the field renormalization even for translation invariant initial data if one considers 
taking kinetic scaling limits of all field observables. We inspect the time correlation of the “bare” 
^-fields, i.e., E['i/:o(k', —l)^t(k, 1)] assuming spatially homogeneous, gauge invariant initial data. 
Then there exists a function 5't(k) such that E['i/)o(k', —l)7/>t(k, 1)] = (I(k' -f k)5't(k). By using 
(5.7) in (4.6), we find the following evolution equation for 'I't: 

(I(k' -f k)^t(k) = «:(V^o(k',I)) = <5(k' + k)4'o(k) - W(k' + k) / dsa;,^(k)4',(k) 

-iX [ ds [ dk2dk3dk4(5(k - ki - k 2 - k 3 )E[:'^o(k',-1): :V't(ki,-l)^t(k 2 , l)^t(k 3 ,1):]. 
Jo 

( 6 . 1 ) 


The linear equation associated to (6.1) thus has the form 


ftik) = /o(k) - i fdsu;^ik)fs{k) , (6.2) 

which is solved by /t(k) = [/t(k)/o(k) where Ut(k) = exp(—i f* dsuj^(k)). We recall that (k) = 
w(k) -I- Ai?s, and thus at a kinetic time scale, with t = tX~^, we have /pdsw^(k) = TA“^u;(k) -|- 
0(A”^). Therefore, Ut has unbounded oscillations in the kinetic scaling limit. Also, we find that 
even though the effect of the first order term proportional to Rg is subdominant, it is still rapidly 
oscillating on the kinetic time-scale and should not be “expanded” in any perturbative treatment 
of the problem. 

We could solve the problem with these unbounded oscillations by considering instead of the 
renormalized field a* = In fact, by the results of Section 5.1.2 and using oq = we find 

that 


<I(k -b k')A,^(k) = E[ao(k', -l)at(k, 1)] = [/t(k)E V'o(k', -l)^i(k, 1) 


S{k + k')Ut{k)^t{k). 

(6.3) 
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We have argued in Section 5.1.2 that the kinetic scaling limit of exists. Then 4'T-;)j-2(k) cannot 
have a convergent limit as A -5- 0; instead, it has fast oscillations proportional to U~^_ 2 {'k). Let us 
also once more stress that the “zeroth order renormalization”, i.e., countering the free evolution, 
does not remove all of the unbounded oscillations but still leaves those resulting from the Rt term. 

However, the above renormalization procedure cannot be straightforwardly extended to more 
complicated cases. Consider next the DNLS model with inhomogeneous initial data. As in section 
5.1.1, our goal is to find the right observable which satisfies the Boltzmann equation in the kinetic 
scaling limit. The evolution equation for the bare field ip reads 

9t^t(k,CT) = -iCTw(k)'^t(k,CT) 

- 2iCTA [ (ikidk2dk35(k - ki - k2 - k3)K(^t(ki,-l),^4(k3, l))^t(k2,CT) 

-iaX [ dkidk2(ik3(5(k - ki - k2 - k3) :^t(ki,-l)^t(k2,CT)^t(k3,1): . (6.4) 

By following the same strategy as for the homogeneous case, let he Ut, t > 0, denote the family 
of linear operators which solves the linear part of (6.4): we suppose that it solves the operator 
equation dtUt = —iHtUt where (iLt/)(k, a) = (Ta;(k)/(k, a) + 2crA dk'Ktpk — k')/(k', a). The 
time-dependent convolution kernel A't(k) should be equal to fjj, d'kiK{ipt(ki, — 1), iptpk - ki, 1)). 
This can be done either by first solving the implicit equation for the above integral over the 
cumulant, or by leaving Kt arbitrary and fixing it by some minimization procedure at the end. 
If such a family Ut can be found, we may define as before at(k, cr) = ({7(“^'0t)(k, ct) and, since 
dtUpp^ = —UU^{dtUt)UU^ = iUpp^Ht, it then satisfies an evolution equation 

a(ht(k, a) = -iaX [ d^k [ dk” [ d^k'Sik” - k[ - k'^ - kf,) 

Jjd J{jd)3 

X {(t(k,k", cr)Mt(k3,ki, -l)Mt(k2,k2,CT)Mt(k3,k3,1) :at(ki, -l)at(k2, cr)at(k3,1):, (6.5) 

where Ut{k, k', a) and Ut(k, k', a) denote the formal integral kernels of the operators Ut and U^^, 
respectively. (Note that the operators Ut are diagonal in cr but not any more in k.) 

Apart from some special cases it seems difficult to gain sufficient control over the operators 
Ut to consider taking a kinetic limit using the observables a*, unlike with the explicit phases 
factors which appeared in the spatially homogeneous case. Even though Ut approach the same 
multiplication operator as before when A —>■ 0 for a fixed t, it is not clear that the corrections do 
not contribute in the limit, since we need to consider t = 0{X~^). Thus, although the cumulant 
expansion of dj-fields is simpler than that of f/'t-fields, to control the kinetic scaling limit looks 
intractable. Hence, new approaches for the study of the kinetic time scales are called for. 


6.2 Kinetic theory beyond kinetic time-scales? 


In this section we propose a new approach to the problem when a renormalization scheme with 
a convergent kinetic scaling limit cannot be found or controlled. The approach does not require 
taking A ^ 0, and, if successful, it may also yield estimates which are valid beyond the standard 
kinetic time scales which was 0{X~^) in the above DNLS case. 

The main idea can be summarized in the following simple observation. Suppose ft is a solution 
to the equation 

ft = fe + Rt+ f dsFs^t[/s] (6.6) 

Jo 

where Rt = 0{e) uniformly in t and Fg^t, 0 < s < t, is an explicit, but possibly nonlinear 
functional of /g. Suppose furthermore that there is another, “simpler”, functional such that 
Fg^t = 4’s,t + 0{e{l -f |t — s|)~^) with p > 1; by simpler we mean that the evolution problem 


(ft = St + / ds$s,t[v3g], 

Jo 


(6.7) 
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for any bounded “source term” St, is easier to study than (6.6). Under these assumptions, any 
solution to (6.6) satisfies 

ft = fo + Pt + f dsd>s_t[/s], (6-8) 

Jo 

where pt := Rt + fgds (Us_t[/s] — $s,t[/s]) is 0(e) uniformly in t. Therefore, if we could prove that 
(6.7) is stable under perturbations of the source term St, we may conclude that the solution (ft to 


Pt = fo+ / ds$s.tbs] 
Jo 


(6.9) 


then approximates the “true” solution ft with an error which is 0(e) uniformly in time. 

To have a concrete example, consider again the DNLS with gauge invariant initial data. We 
sketch below two conjectures about this system. The details of the conjectures should not be 
taken too seriously: they should be considered more as examples of what could happen in general 
rather than as specific conjectures about the behavior of the DNLS system. For this reason, the 
discussion will be kept on a very loose level; in particular, we do not wish to make any specific 
claims about what kind of metrics should be used for studying the uniform boundedness in time. 

Let us begin with the case of inhomogeneous, gauge invariant initial data. We move to slowly 
varying fields by cancelling the free evolution term. This renormalization leads to equations 
which are almost identical to those in the homogeneous case. Namely, for the field &t(k, ct) := 
ei‘^“(k)*'i/i((k, cr) we have 


dMKa) = -2iaA / dkidk 2 dk 3 ( 5 (k - ki - k 2 - k 3 )e*‘("‘"+'"i--‘"^-‘" 3 )Bt(ki, kg) : 6 i(k 2 , a): 
J{-^'tY 

- iaA f dkidk 2 dk 3 (j(k - ki - k 2 - k 3 )eid--+-i—2-^3) -l)&t(k 2 , a)bt(k3, 1): 

d(T‘^)3 

( 6 . 10 ) 

where i?t(k', k) = K( 6 t(k', —1), 6 t(k, 1)). By the gauge invariance, K(&((k', ct'), 6 t(k, a)) = 0, unless 
ct' + cr = 0. Hence, there is only one other nonzero second order cumulant, K(bt(k', 1), 6 t(k, —1)) = 
Bt{k,k') = Bt(-k',-k)*. 

To study the fourth order cumulants, it suffices to concentrate on the function 
Dt(k) := K(&t(ki, -1), &t(k 2 , -1), 6 t(k 3 ,1), 6t(k4, l)),k€ (jy . 


Then the first two equations in the cumulant hierarchy are equivalent to 

dtBt(k',k) = -2iA [ dkidk 2 dk 36 (k - ki - k2 - k3)e'*(‘"+‘"i-‘"^-‘"^)Bt(ki,k3)Ht(k2,k') 

+ 2iA f dkidk2dk3S{k' - ki - k2 - k3)e'*(-‘"'+‘"i+‘"^-‘"^)Ht(ki, k3)Ht(k, k2) 

J{'j't)3 

- i\ [ dkidk2dk36ik - ki - k2 - k3)e‘‘(“+‘"i-‘""-‘"3) A(k', ki, k2, ka) 

+ i\ [ dkidk2dk35(k'-ki-k2-k3)e'*(-‘"'+‘"i+‘"=-“^)A(ki,k2,k3,k), (6.11) 

J{1<i)3 

4 /. 

dtDt(k) = -i2X^ae / dk'3dk^dk(,(j(k£ - k^ - k^ - 

X Ht(k 3 ,k 3 )A(“replace ke by k^”) 

^ r , , , 

- iA V / dk;dk^dk(,(j(k£ - k; - k^ - 

X (“(1 X y + (5 X BtDt) + (6 X BtBtBt )”), (6.12) 
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where in the second formula, a := (—1, —1,1,1) and on the last line we have applied (5.22) and 
merely denoted how many nonzero terms each type of partition can have. Let us point out that 
the first two terms in (6.11) cancel each other out if the state is spatially homogeneous, since then 
ili(k', k) oc d(k' + k). For an inhomogeneous state, however, the cancellation need not be exact, 
and since this term is then not O(A^), it will likely prevent taking of the kinetic scaling limit of 
Bt directly. 

If we now integrate the B equation as in (4.12), we find that 


Bt(k', k) = Bo(k', k) - iA / ds / dkidk 2 dk 3 d(k - ki - ka - 


X (2i?s(ki, k3)i?s(k2, k') + iAo(k', ki, ka, k3)) 


+ iA / ds / dkidk 2 dk 3 d(k'- ki - ka - k 3 )e'®(-‘^'+“i+‘^"-“=^) 
Jo JiJd-Y 

X (2Bs(ki,k3)Bs(k, ka) + iAo(ki, ka, k3, k)) 

+ \^ f ds' f ds ^ (•••)• 


(6.13) 


In the final sum, each term depends on s only via the oscillatory phases. These can be collected 
together and they have a structure 


^is{uo{u:o—‘^2}+^^i—‘^3)Qis'{aeu>e+u:[ — aeu:'2—uj'^) 


(6.14) 


where (CTo,ko) is equal to (l,k), if the term arises from the second last term in (6.11) and it 
is equal (— l,k') if it arises from the last term. The pair comes from arguments of the 

corresponding “Ht-term” and thus depends also on this choice. Therefore, the s-integral over the 
oscillatory phase can be computed explicitly and to each of the terms in the sum it will produce 
a factor 




nt — s' 


dr e 


ir{ao{c^o—^2)+oJl—‘^3) 


(6.15) 


It is difficult to go further in the analysis of the oscillatory phases without resorting to graph 
theory, and we will not pursue it here. However, already the simple example given in Appendix 
C shows that, if the state is fi-clustering, the oscillations may result in time-integrals which are 
absolutely convergent over [0,cx)). For instance, this explicit example implies that if the initial 
state is homogeneous and £i-clustering, then the two terms depending on Dq in (6.13) are 0(A) 
uniformly in time. If we assume that similar bounds are valid for inhomogeneous states and every 
term containing kq, the two-component field ft = {Bt, Dt) behaves as the model considered in the 
beginning of the section. 

Therefore, assuming uniform boundedness of Kg allows using the principles described in the 
beginning of this section and results in a conjecture about the evolution of the cumulants. Suppose 
that there is a metric for the cumulants, similar or given by fi-clustering estimates, such that the 
following results hold for sufficiently regular dispersion relations w and initial data: 

1. Suppose that the contribution from KQ{t) is uniformly bounded in time, with a bound 
0{X^~^) where 0 < g < 1. 

2. Consider the evolution equation obtained for {Bt,Dt) by the standard “closure relation”, 
i.e., by setting Kg —> 0 in (6.12). Assume that this equation is stable under all uniformly 
bounded time-dependent perturbations of the source term. 


If both of the above hold, then the solutions to the closure evolution equation for ft remain 
close to {Bt,Dt) uniformly in t. In particular, the difference in the first (covariance) component is 
always 0(A^“*). This would validate the closure equations as good approximations even to f —5- oo 
asymptotic behavior of the covariance function for all sufficiently small A. Note that no scaling 
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limit needs to be taken; in particular, it is not claimed that the kinetic scaling limit of Bt or Dt 
would exist. 

Another application can be obtained for a one-component case with ft = Wt as follows: 
Consider the spatially homogeneous case and the exact evolution equation for obtained from 
(5.18) by “dividing out” (5(k'-|-k) from both sides. Take Rt to include all terms which contain either 
K 4 or Kg- Then the remaining pairing terms yield an explicit definition for Fg^t such that (6.6) 
holds. (In fact, then Fs^t\W] is equal to the right hand side of (5.26).) Next choose equal to 

the Boltzmann collision operator A^C with C defined in (5.14). If Wt comes from an £i-clustering 
state and the free evolution is sufficiently dispersive, then Fs^t\Ws\ = C[lTs] -I- 0(A^(1 -f |t — s|)“^), 
with p > 1. (For instance, the estimates given in [7], in Proposition 7.4 and in the Appendix, 
prove the bound with p = 3d/7 — 1 for a nearest neighbor dispersion relation—the computation is 
essentially the same as in Appendix C below. Using the notations defined in (C.3), the result also 
allows to quantify the dependence on the £i-clustering assumption: the bound is proportional to 
||'« 2 (s)||f A^(l + \t — s|)“^. Hence, then p > 1 at least if d > 5. However, more careful estimates or 
the addition of next to nearest neighbor hopping could improve the bound.) 

Whenever this is the case, we obtain a second conjecture about the homogeneous DNLS equa¬ 
tion. Suppose that there is a metric for the cumulants such that the following results hold for 
sufficiently regular to and initial data: 

1. Suppose that the cumulants remain uniformly bounded in this metric, with an upper bound 
which implies that ||k 2 ||i = 0(1) and that all higher order cumulants have -clustering 
norm which is 0(A“*) for some 0 < g < 2. 

2. Assume that the corresponding Boltzmann-Peierls equation is stable under all uniformly 
bounded time-dependent perturbations of the source term. 

If both of the above hold, then the solutions HA- to the Boltzmann-Peierls equation with initial data 
Wq are 0(A^“'^) close to W ^^_2 uniformly in r. In particular, any stationary limit limj^oo Wt 
can differ from the limit of the solution to the Boltzmann-Peierls equation only by 0(A^“®). 

The main benefit from using the Boltzmann-Peierls equation instead of the closure hierarchy 
concerns the second assumption: the homogeneous Boltzmann-Peierls equations enjoy many sim¬ 
plifying properties and a priori estimates, see for instance [9]. For example, they typically have an 
entropy functional and an associated “H-theorem” which allow to classify all stationary solutions 
to the equation. There are also many techniques developed to control the convergence towards 
the stationary solution. 

As a final example, let us remark that even when total uniformity in time cannot be achieved, 
it might be possible to go beyond the kinetic time-scales using the above methods. Consider 
Fs,t\Wt] = C\Wt] 0(A^(1 -f |t — s|)~^) for some 0 < p < 1. Then the correction is not in- 
tegrable and the perturbation pt to the source term is 0(A^t^“^). This remains 0(A®) for all 
t = Hence, for instance, the earlier nearest neighbor estimate with d = 3 would 

imply that the corrections to the Boltzmann-Peierls equation remain small for t = 
that is, even for times much longer than the ones implied by the kinetic scaling limit. 

A Combinatorial definition of cumulants 

It is shown in [8] that cumulants are connected to moments via a formula which is very similar to 
the definition we used here for the Wick polynomials: if / 0, for any a; S / we have 

E[/] = ■ (A.i) 

E-.x^Edl 

(The formula follows straightforwardly from the identity dxGm = Gjndxdc) In fact, this formula 
allows a definition of cumulants which does not rely on differentiation or on the existence of 
exponential moments. Namely, ii Iq € is such that E[|2/'^|] < oo for all / C Iq, then to each I C 
lo, I 7 ^ 0, we can associate a number K[yi] by requiring that K[yi] = ^[y^]—'^E-xeEci ^[y^'^^]i^[yE] 
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with X = The definition is used inductively in |/| > 1 and it has a unique solution. (Note 

that the empty cumulant ^[ 2 / 0 ] never appears in the moments-to-cumulants formula, and for our 
purposes it can be left undefined. To be consistent with the derivatives of the generating function, 
we may for instance set ^[ 2 / 0 ] := 0.) Therefore, these numbers have to coincide with the standard 
cumulants in the case when exponential moments exist and hence (A.l) holds. 

The following known properties of cumulants can then be derived directly from the above 
definition using induction in |/| and applying techniques similar to what we have used for Wick 
polynomials in Section 3.1: 

1. The cumulants are multilinear, in the same manner as was stated for Wick polynomials in 
Proposition 3.7. 

2. The moments-to-cumulants expansion (3.1) holds. 

3. The cumulants are permutation invariant: if I' is a permutation of I, then «:[ 2 //'] = n[yi\. 

4. If joint exponential moments exist, then K[yi\ = 9{2/c(0) with ( 2 c(A) := lnE[e^'*^]. 

However, let us skip the proofs here. In the text, we assume these results to be known and refer 
to the references for details of their proofs. 

B Classical particle system with random initial data 

Consider the evolution of N classical particles interacting via a polynomial interaction potential, 
with the initial data given by some random probability measure. We show here how it can be 
recast in the form of the evolution equation discussed in Section 4. 

We consider the random variables yj, indexed by J = {(z,n)}i_n, where n is one of the N 
different particle labels and f = 1,2 differentiates between the particle position and momentum: 
we dehne 2/(i,ra)(0 := and 2 /( 2 ,n)(^) := Pn{t)- If all particles move in OR, have the same mass, 
and have only pair interactions via the potential V{q) := 'Yhn' n-n'^n Ann' ^(^n — <?n')“, a > 2 even 
and A„'„ = A„„', then we have dtqn{t) = Pn{t) and 

dtPn{t) = - ^ Kn'{qn{t) - g„'(f))““^. (B.I) 

n'^n 

Here (g„ - and if we define In,n',k as a sequence of 

length a — 1 containing first k repetitions of (1,n) and then a — l — k repetitions of (I, n'), then by 
(3.6) we have qn{tYq-n'^ £(2/(0^"’"'’'“^^)- Define thus as the collection 

=^( 2 .n) all such sequences U which contain ki repetitions of (I,n) followed by ^2 repetitions of 
(1, n') where n' yf n, fci, ^2 > 0, and fci -I- ^2 < a — I. Set also for each U € -*^( 2 ,n) 

A^(^.n)W:=Ann' ^ E E(2/1 (D = t/). (B.2) 

Therefore, (4.3) holds for all j G J after we also define J^(i,n) := {0, ((2,n))} and set ^^(t) := 
Hy( 2 ,n){t)) and := 1 if [7 = ((2,n)). 

C Estimation of the first order non-pairing contributions 
to (5.18) 

In this appendix, we show how to estimate the first order non-pairing contraction terms in (5.18). 
To this end we need to make an assumption on the dispersion relation w(k). Let us consider the 


27 



so called free propagator 


Pt{^)= [ . (C.l) 

Jjii 

As in the assumption “(DR2)” in [7], we now suppose that there are C,S > 0 such that for all 
t ^ IR, 


|b,||i= ^ b*(x)|3<C(l + t2)-(i+^)/^ 


(C.2) 


Furthermore, we assume also the already mentioned t^i-clustering property (see section 4) which 
we slightly rephrase as follows for each cumulant of order n: we require that 


||k„||i := sup ^ l(xi = 0)|K(V’(xi,cri),...,V’(x„,CT„))| < oo. 


o-e{±i} 


(C.3) 




We recall that the physical meaning of this condition is that the cumulants decay fast enough in 
space so that they are summable, once the translational invariance is taken into account. 

We recall that the first order non-pairing contributions in (5.18) are 


iAu f ds f dkidk 2 dk 3 ( 5 (k - ki - k2 - 
X K[a(ki, -1); a(k2, ct); a(k3,1); a(k', a')] 


(C.4) 


and a term which is obtained from (C.4) by swapping (k,cr) -O- (k'jCr'). As stated in (4.1), by 
translation invariance we have 

K[a(ki, - 1 ); a(k2, cr); a(k3, 1 ); a(k', cr')] = 5 (ki -|-k2 -b ks -b k')F(ki, k2, ks, k', ct, ct') (C. 5 ) 

where F'(xi, X 2 , X 3 , X 4 , cr, cr') = l(xi = 0)K[a(xi, -l);a(x2,CT); a(x 3 ,1); a(x 4 , ct')]. Clearly, ||F||i < 
||k 4 ||i < 00 by the assumed £i-clustering. 

Therefore, the term in (C.4) is bounded by 


A 


f ds [ dkidk2dk3(5(k - ki - k 2 - k3)e'"('""'+“i-'""'"-"'3) 

Jo J(j<i)3 

X ( 5 (ki -bk2 -bk3 -bk')F(ki,k2,k3,k',cr, cr') 


< 


f dkidk2 e‘*(‘"i-""'=-“^)i?(ki,k2,k3,k',a,CT') 


d(T‘i)2 

k3=k-ki-k2 


< A(5(k -b k') / ds E |F(xi,X2,X3,X4,ct, a'] 

0 Xi,X2,X3,X4 

Ee“i'"‘^'(^”x^)p-^(yK.(y - X2)p.(y - X3) 

y 

< A(5(k + k')||K4||i f ds Ijpsllg < AC'5(k + k')||K4||i f ds(l + 

Jo Jo 

< AC'd(k + k')||«; 4 ||i 


(C.6) 


where C is a constant which depends only on C and 5. We have used the inverse Fourier transform 
of F and (C.l) in the second inequality and Holder’s inequality in the third one. Since the bound 
in invariant under the swap (k, cr) o (k',CT'), it bounds also the second non-pairing contribution 
in (5.18). Therefore, we see that the first order contributions are 0(A) uniformly in t. 
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